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ABSTRACT 

We have observed a small sample of hot post-AGB stars with the InfraRed Array Camera (IRAC) 
and the InfraRed Spectrograph (IRS) on-board the Spitzer Space Telescope. The stars were selected 
from the literature on the basis of their far-Infrared excess (i.e., post-AGB candidates) and B spectral 
type (i.e., close to the ionization of the envelope). The combination of our IRAC observations with 
2MASS and IRAS catalog data, along with previous radio observations in the cm range (where avail- 
able) allowed us to model the SEDs of our targets and find that in almost all of them at least two 
shells of dust at different temperatures must be present, the hot dust component ranging up to 10^ K. 
In several targets grains larger than 1 /xm are needed to match the far-IR data points. In particular, 
in IRAS 17423-1755 grains up to 100 /xm must be introduced to match the emission in the mm range. 

We obtained IRS spectra to identify the chemistry of the envelopes and found that more than ^/a of 
the sources in our sample have mixed chemistry, showing both mid-IR bands attributed to Polycyclic 
Aromatic Hydrocarbons (PAH) and silicate features. The analysis of the PAH features indicates that 
these molecules are located in the outflows, far away from the central stars. We consider the larger 
than expected percentage of mixed-chemistry targets as a selection bias towards stars with a disk or 
torus around them. Our results strengthen the current picture of mixed chemistry being due to the 
spatial segregation of different dust populations in the envelopes. 

Subject headings: Infrared: general, Stars: AGB and post-AGB, Planetary Nebulae: general 



1. INTRODUCTION 

Planetary Nebulae (PNe) evolve from intermediate 
mass stars. During their evolution they go through the 
Asymptotic Giant Branch (AGB) and then into the PN 
phase, ending their days as white dwarfs. The forma- 
tion and early evolution of PNe are a poorly understood 
phase of stellar evolution. In particular, it is not clear yet 
how the spherically symmetric AGB circumstellar shells 
transform themselves into the non spherical symmetries 
observed in the envelopes of evolved PNe. During their 
AGB phase stars lose most of their initial mass, their 
mass loss rates ranging from 10^^ up to lO"** M0 yr~^. 
The star develops a wind that carries away from its pho- 
tosphere the enriched material produced by nucleosyn- 
thesis and dredged up from the core to the surface. This 
ejected material eventually cools down and creates a cir- 
cumstellar envelope of dust and gas, whose chemistry is 
determined by the C/0 ratio in the stellar outer regions 
([Hiriart & Kwan 2000). When the star has evolved past 
the tip of the AGB, the dust and gas circumstellar enve- 
lope is typically so thick that the star cannot be detected 
at optical wavelengths: the expansion and subsequent 
dilution of the envelope leads to the optical detection of 
post-AGB stars. The end of the post-AGB evolution can 
be identified as the onset of the ionization in the enve- 
lope, when the central star reaches temperatures around 
20-30x10'^ K. The ionization is in fact considered as the 
beginning of the PN phase. 

Studies of the infrared (IR) properties of AGB and 



post-AGB stars have determined that in 0-rich envelopes 
silicate features at 10 and 18 /Ltm are first seen in emission 
during the AGB, then in absorption when the envelope 
becomes optically thick, and then again in emission, after 
the en velope becomes optically thin to the m id-IR radi- 
ation l|Garcia-Lario fc Perea CalderonI 120031 ). A similar 
evolutionary scheme is observed in C-rich envelopes, al- 
though with different spectral features. Low-mass 0-rich 
stars are not expected to develop an optically thick enve- 
lope, therefore the silicate features a re never seen in ab- 
sorpt ion in these stars (Garcia-Lario fc Perea CalderonI 
[200l . 

A small fraction of post-AGB stars and PNe show both 
0-rich and C-rich features and are therefore classified 
as mixed- chemistry objects. The origin of their mixed 
chemistry is still uncertain. For example, in the Red 
Rectangle this has been explained as due to the presence 
of a circumbinary disk that was formed when the star was 
still 0-rich and then served as a silicate reservoir, while 
the C-rich features originate in the dredged-up outflow 
(jWaters et al.lll998r ). A different scenario has been pro- 
posed for high mass stars, in which hot bottom burning 
would convert ^^G into ^"^ N and change a C-rich envelope 
back into an 0-rich one (jJustanott et al.l[l992[ ). 

Many observational programs have been carried out 
to find new planetary and proto-planetary nebulae 
among unidentified IRAS sources with far inf rared col- 
ors similar to those of known PN e (|Pottasc h~et al. 198^; 
iParthasarathv k, PottaschI I1989D . with the final aim of 
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understanding PN formation through the discovery and 
analysis of new Transition Objects (TO). In these sources 
the physical processes associated with PN formation, 
such as dynamical shaping, are still occurring and, as 
the number of known TO is extremely small, the identi- 
fication and study of new samples is relevant for further 
testing current models of stellar evolution. 

To detect new TO, we have selected a sample of stars 
classified in the literature as post-AGB on the basis of 
their IRAS colors and optical spectra. We have selected 
our stars from li sts of post-AGB candidates available 
in th e literature (jParthsarathv et al.ll2000t ISuarez et al.l 
120061) ■ with the constraint of a B-type spectral classifi- 
cation and far-IR excess. The B spectral type implies 
the central star can be hot enough to have started to 
ionize the circumstellar envelope. This maximizes the 
chance that the objects are right in the transition from 
post-AGB to PN, although misclassifications cannot be 
ruled out. Radio observations of this sample were carried 
out to check for the presence of free-free emission from 
ionized shells . Thes e obs ervations have b e en repo rted in 
lUmana et all (|2OO40 and lCerrigone et al.l (|2008al ). 

To further characterize our sample, we performed ob- 
servations with the Spitzer Space Telescope. Mid-IR ob- 
servations probe the extended envelopes that surround 
post-AGB stars and young PNe, which are in most cases 
only partly ionized. The main goals of such observations 
are to assess the morphology of the neutral gas and dust 
components and to characterize their mineralogy. Com- 
paring the results of these observations to the properties 
of the ionized g derived from previous radio obser- 
vations, could in principle provide important clues on the 
shaping of the nebulae. 

2. OBSERVATIONS 

As a follow-up to our radio detections, we observed 26 
objects with the Spitzer Space Telescope (Spitzer). In- 
frared imaging w as performed with the InfraRed Array 
Camera (IRAC) ()Fazio et al.ll2004[ ) at 3.6, 4.5, 5.8, and 
8.0 /im to search for faint extended emission. The ob- 
servations were carried out in several Spitzer campaigns 
in 2006 and 2007 (Program ID: 30036). The IRAC 
observations were performed in High Dynamic Range 
mode. For each target we obtained 8 dithered frames, 
for a total on source time of about 96 s per IRAC chan- 
nel. Basic Calibrated Data (BCD) were retrieved from 
the Spitzer archive (pipeline versions: S14.0.0, S15.0.5, 
and SI 6. 1.0), cleaned to correct such artifacts as mux- 
bleeding and banding and then coadded using IRACproc 
(jSchuster et al.ll2006f) . The photometry was performed 
by PSF subtraction in the final mosaics. This tech- 
nique has been specifically developed to p erform photom- 
etry even in heavily saturated objects ([Marengo et al.l 
I2007| ). The method derives Ve ga magnitudes of saturated 
stars by fitting the low level features of the saturated 
PSF (diffraction spikes and PSF wings) with a model 
of the IRAC PSF obt ained from a sample of bright stars 
(jMarengo et al.|[200^ . The PSF is directly normahzed to 
the IRAC observation of Vega, therefore, by fitting it to 
the observed sources, magnitudes can be calculated with 
an accuracy within 1- 3% independently of the standard 
IRAC flux calibration (jSchuster et al.ll2006D . 

Observations w ith the InfraRed Spectrograph (IRS; 
iHouck et al.ll200^ were performed within the same pro- 



gram in 2006 and 2007. The observations were carried 
out in low spectral resolution mode (R~60-120), cover- 
ing the 5-38 /xm range, to assess the chemistry of the 
envelopes. Each target was observed in all the available 
low resolution modules (SL and LL) with a ramp time of 
6 s per cycle. Two cycles per module were typically per- 
formed. Those targets with IRAS 12 /j,m flux below 1 Jy 
were observed for 3 cycles in the 5.2-14.5 /im modules, 
to keep a high enough S/N (typically ~50). The IRS 
pipeline version S17.2.0 (partial) BCDs were retrieved 
from the archive. The sky subtraction was performed 
taking advantage of the way the data are acquired by 
default. Each module has two sub-modules (SL2 and 
SLl, LL2 and LLl) observing at the same time, but not 
in the same direction. Since the two sub-modules do not 
point the same target at the same time, when one is on- 
target, the other is performing a sky observation. The 
data taken by SLl (LLl) when SL2 (LL2) is on target 
can then be used as a sky frame for a subsequent SLl 
(LLl) pointing and vice versa. For each BCD, we have 
averaged all the off-target observations and subtracted 
this from the single on-target BCDs. The spectrum was 
extracted in every BCD with SPICE using the optimal 
point source extraction option. For each target, all its 
extracted spectra we re then averaged an d combined into 
one using SMART^ ()Higdon et al.ll2004[ ). 

3. RESULTS 
3.1. IRAC photometry 

IRAC has proved to be a very sensitive instrument, re- 
vealing previously unseen structures in the outer shells 
of expanded PNe ( Hora et al, 2004 ). Since the nebulae 
in our sample are typically very compact and IR-bright, 
which is normal for their evolutionary status, their core 
emission dominates the inner few arcsec near the central 
source. This unfortunately usually prevents us from de- 
tecting any weak extended emission near the core, but, 
by PSF subtraction, does not prevent us from measur- 
ing the fluxes of our sources. Given the size of the PSF, 
an upper limit of '^5" to the sizes of the nebulae can be 
estimated (Marengo, private communication). 

In Table [1] we list the photometric results ob- 
tained. For a comparison of our results to p revious 
observations in the liter a ture (i.e., iHora et al.l ((2008,) : 
IPhillips fc Ramos Larioj (|200l) ). we show in Figure [T] 
four color-color diagrams b ased on IRAC data. As al- 
ready shown in iHora et al.l (2008), PNe do not seem to 
follow any speciflc track on the IRAC [3.6] — [4.5] ver- 
sus [5.8] — [8.0] diagram. We find our targets to be gen- 
erally located in an area ranging from to 3 in [5.8] — [8.0] 
and from to 1 in [3.6] — [4.5]. This region is included in 
the area where PNe are found in the two works mentioned 
above. A comparison of the [3.6 | — [4.5] versus [4.5 ] — [8.0] 
diagram to the analogous one in] Hora et all (|2008[ ) shows 
that the stars in our sample have on average lower values 
of [3.6] — [4.5], which may be explained if PAH emission 
features are present (see S I3.2p : the spread in [4.5] — [8.0] 
indicates that different dust temperatures can be found 
in these objects, with the radio-detected ones (i.e., hotter 
central stars) being redder in [4.5] — [8.0]. 

^ SMART was developed by the IRS Team at Cornell University 
and is available through the Spitzer Science Center at Caltech. 
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TABLE 1 

IRAC FLUXES. 



Target 


F3.6 


F4.5 


F5.8 


Fg.o 


F24 


IRAS ID 


mJy 


mJy 


mJy 


mJy 


Jy 


(JiUU5+7yiU 


175. D ± 9.1 


1 on 1 n 
I59.Z ± 9.5 


460. U ± 24. U 


onn At 1 -y n A A 

20U4.1 ± 1U4.4 


18.1 ± 0.9 


Uooob+ibzo 


vc^n _L on ^ 
(0\).Z ± oy.o 


595.0 ± 4d.o 


1U40.0 ± 04.0 


1 1 n n _L cii n 
121U.U ± 00. U 


2.8+ 0.2 


nn A '7r\ A£ii '7 


rri _L Q '7 
IV.Z ± 0. ( 


on n _L A 7 
59.9 zt 4. ( 


1 n 1 7 _L inn 
191./ ± 10. U 


12oo.o ± d4.o 


3.5 + 0.2 


17U74-is45 


17.6 ± 0.9 


"ion 1 n £? 

12.0 ± 0.6 


10.5 ± 0.5 


47.5 di 2.5 


in'? 1 n 

1U.7 ± 0.6 


172Uo-i5ci4 


18.7 ± 1.0 


100 1 n T 

12.8 ± 0.7 


11.5 ± 0.6 


24.7 di 1.3 


8.9 ± 0.5 


17viD4-i2ci» 


lu.o ± 0.6 


8.2 ± 0.4 


7.7 ± 0.4 


9.9 di 0.5 


1.1 ± 0.1 


17jOl-iDlD 


4.3 ± 0.2 


3.6 ± 0.2 


2.9 ± 0.1 


33.1 di 1.7 


4.4 ± 0.3 


17423-1755 


3425. D ± 175.4 


4492.5 ± 2v>o.5 


5227. ± 272.5 


c /I H 1 or^o cr 

5o44.2 ± 275.5 


r\ A 1 TO 

24.2 di 1.2 


1 /o4z-UdUo 


nOC Q _L /ICQ 


1U5/.1 ± 55. U 


1U45.5 ± 54.5 


line; 7 _L Ti'7 a 
11U5. ( ± 0( .D 


5.3 di 0.3 


loU4U-145 / 


Q 1 1 _L 1 Q 


01 1 A _L 1 1 n 
211.4 lb 11. U 


onn 1 _L inn 
2U9.1 ± 1U.9 


1 /1 c; _L 7 A 
142.5 ih /.4 


O.lldiO.ll 


1oUd2+241U 


7.4 di 0.4 


9.0 di 0.5 


13.5 di 0.7 


557.7 ± 29.1 


-\ c '~> 1 n n 

16.3 ± 0.9 


1oU7U-z^s4d 


1 '7 1 no 

lo7.o ± 9.5 


1 1 n 1 r" 

119.8 ± D.2 


n 1 1 /I 

82.1 ± 4.3 


55.8 ± 2.9 


1.0 ± 0.1 


lOOD /-1200 


DOOZ.D ± o4U.O 


rrinn n _l_ qi i o 
0990. U ± oll.o 


c^'zc^n n _L onn o 
/OO.U ± 299.8 


/ 125.6 ± ( 1.0 




i8vi7i-oi5y 


7.0 dl 0.4 


6.0 ± 0.3 


7.7 ± 0.4 


22.1 ± 1.2 


5.9 ± 0.3 


loo / y-1 / U ( 


41. it 2.Z 


00. z it z.U 


oz.i it 4.0 


Qon 7 -L 1 7 

OZU. ( it 10. ( 


on n 1 1 n 

zy.y it i.u 


18442-1144 


16.5 ± 0.9 


18.0 ± 0.9 


71.9 ± 3.7 


320.7 ± 16.7 


12.1 ± 0.6 


19200+3457 


18.7 ± 1.0 


13.8 ± 0.7 


38.3 ± 2.0 


142.5 ± 7.4 


1.6 ± 0.1 


19306+1407 


40.1 ± 2.1 


30.0 ± 1.6 


88.5 ± 4.6 


583.0 ± 30.4 


47.5+ 2.9 


19336-0400 


9.4 ± 0.5 


9.0 ± 0.5 


10.5 ± 0.5 


53.4 ± 2.8 


6.7 ± 0.4 


19399+2312 


280.9 ± 14.6 


179.7 ± 9.4 


164.3 ± 8.6 


142.5 ± 7.4 


0.12+0.11 


19590-1249 


14.0 ± 0.7 


11.2 ± 0.6 


11.5 ± 0.6 


25.7 ± 1.3 


7.9 + 0.4 


20462+3416 


24.4 ± 1.3 


16.3 ± 0.9 


17.7 ± 0.9 


27.9 ± 1.5 


11.1 + 0.6 


20572+4919 


936.3 ± 48.8 


898.5 ± 46.8 


884.6 ± 46.1 


1733.2 ± 90.3 


8.5 + 0.4 


21289+5815 


61.1 ± 3.2 


64.2 ± 3.3 


92.0 ± 4.8 


246.7 ± 12.9 


1.1 + 0.1 


22023+5249 


18.7 ± 1.0 


16.3 ± 0.9 


37.7 ± 2.0 


229.0 ± 11.9 


19.0 + 1.1 


22495+5134 


10.4 ± 0.5 


9.5 ± 0.5 


12.1 ± 0.6 


80.2 ± 4.2 


9.7 + 0.5 



Note. — The values in columns 2-5 are calculated by PSF subtraction. The errors include 
the error in the IRAC zero point magnitudes and a 5% error due to the PSF subtraction method. 
The values at 24 /xm listed in column 6 are calculated from IRS spectra. 

No IRS spectrum available to calculate the 24 ^m flux density. 

We also plotted our targets on a H-K versus K-[3.6] 
diagram, which can be compared to the JHK diagram in 
FigurelU A similar trend can be fou nd in both plots. As 
shown in Ivan der Veen et"al] (|1989[) . the JHK dia gram 
allows us to distinguish different types of sources. Many 
of our targets are concentrated in the hot star and neb- 
ula + star regions, as expected for their B spectral type 
classification. Other targets fall within regions of objects 
with a hot star and strong dust emission. The presence of 
a large amount of dust around the central stars is also ev- 
ident in the [3.6] — [8.0] versus [8.0] — [24] diagram. Since 
we do not have MIPS observations of our targets but we 
have 5-38 /im IRS spectra, we calculated the expected 
24 /im MIPS measurement by convolving our IRS spectra 
with the MIPS Spectral Response Function and then in- 
tegrating over the bandpass. Performing this calculation 
on literature targets with both IRS and MIPS data has 
produced results compatible, within errors, with the ob- 
servations, therefore we are confident that our expected 
MIPS values are reasonable approximations to the actual 
values. Interestingly, among the radio detected objects in 
this diagram we can distinguish a group of stars clustered 
around (4,2), namely IRAS 06556+1623 and 17423-1755, 
while the others seem to be found in a stripe roughly 
going from (5,6) to (9,1). After noticing these possible 
different distributions, we checked the positions of the 
"central" targets in the other diagrams, thus finding that 
although the ove rall distribu t ions o f our targets match 
those observed in lHora et all (|2008[ ). IRAS 06556+1623 
and 17423-1755 seem not to follow the trend of the other 
radio detected objects. In [3.6] — [4.5] versus [5.8] — [8.0] 
and [3.6] — [4.5] versus [4.5] — [8.0], these stars are found 



respectively around (0.7,0.7) and (1.3,0.7), while in H- 
K versus K-[3.6] they are the two top- right radio detected 
objects (highest H-K and K-[3.6]). 

3.2. IRS spectra 

The reduced IRS spectra of the objects in our sam- 
ple are shown in Figures [SHZl They are split into three 
groups: C-rich, 0-rich, and mixed chemistry. PNe are 
usually classified as C-rich or 0-rich. This classification 
is linked to the evolution of the central star, because it 
depends on the possibility for the star to go through a 
third dredge-up. It is this event that alters the chemistry 
in the envelope t urning it from 0-rich in to C-rich during 
the AGB phase (jSalaris fc Cassisill2005f ). Because of the 
stability of the CO molecule, if C is less abundant than 
O, all C is trapped in CO and then the envelope shows 
features of molecules containing oxygen (besides CO) and 
vice versa if O is less abundant than C. It is known how- 
ever that a minority of stars shows mixed chemistry, with 
C-bearing and 0-bearing molecules. 

More than ^ of the stars in our sample show both 
Polycyclic Aromatic Hydrocarbons (PAH) and silicate 
features, while a smaller percentage (^^10%) is expected. 
Among these mixed-chemistry objects about 70% are ra- 
dio detected. 

We distinguished 0-rich and C-rich envelopes on the 
presence of amorphous silicate (10-/im bump) or PAH 
features (mainly 6.2 and 7.7 /im). Objects with both a 
10-/im silicate feature and PAH transitions were classified 
as mixed-chemistry. It must be noticed that the SL and 
LL IRS modules overlap around 14-15 /im, therefore this 
region of the spectra is subject to higher noise. Shifts 
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Fig. 1. — Distributions of the sample in four color-color diagrams. The sources are represented as squares if detected at radio wavelengths, 
as diamonds if not. 

All the targets classified as C-rich show 6.2, 7.7, 8.5, 
: — ' — I — ' — ' — ' — I — ' — ' — ' — I — ' — ' — ' — q and 11.2 /zm features, commonly attributed to PAH. The 

shape and strength of the features depend on the exci- 
tation status and size of the molecules. We address the 
properties of these features in the next section. 

We do not detect the so-called 21 feature in any 
of our targets. This feature has in fact been found in 
15 po st-AGB stars with F-G spectral types and in two 
PNe l! Molster fc Watersll2003 l ) . Therefore, besides being 
a rare feature, the chances for its detection are even lower 
in B stars as those in our sample, because of the harder 
radiation field from the central star. 

The sources IRAS 01005-^7910 and 19200-^3457 show 
a strong broad feature around 30 fxin that has been pre- 
viously observed in other post-AGB stars and PNe and 
is attributed to MgS. The peak position and width of 
this feature change as the temperature of the underly- 
ing continuum changes, corresponding to the evolution 
from post-A GB to PN. The pea k is found in the 26- 
35 fim range (jMolster et al.|[200^ . It is likely that such 
emission is present in IRAS 18442-1144 around 34 fj,m, al- 
though weaker and with a smaller width. This target and 
IRAS 06556+1623 show an excess of emission peaking 
around 17.7 /im (Figure[5|) that can b e explained as due 
to vi brationally excited C-H bonds (jMolster fc Wateri 
I2003f) . This emission is also detected in 01005-^7910 and 
19200-^3457 as a very broad feature. IRAS 06556-f 1623 
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Fig. 2. — Distribution of the sample in the JHK color-color di- 
agram. As in Fig. [T] the sources are represented as squares if 
detected at radio wavelengths, as diamonds if not. Radio detected 
targets tend t o concentrate in what has been identified as the hot 
star region in Ivan der Veen eraLl ||T989D . 

between the flux levels are also possible. In a few cases, 
we have scaled the SL spectra to the LL flux level. 
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Fig. 3. — IRS spectra of sources classified as C-rich on the basis of their PAH features. IRAS 1742 3-1755 is listed in this group because, 
though not showing any PAH feature, it shows C2H2 absorption feature l|Gauba fc Parthasarathvl 2004). The whole spectrum of each 
object is shown on the right. On the left we show the 5-14 fim interval of the same spectrum after the subtraction of a locally-defined 
continuum, to enhance the features detected in this wavelength range. As an example, error bars are overplotted for IRAS 19200-1-3457. 
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Fig. 5. — As in Figure |4] 



also shows a weak broad emission underlying the 6.2, 7.7 
and 11.2 /im features, which is due to PAH molecules as 
well. 

In 0-rich targets, as we can see in Figures |4] and [5l the 
10 /im broad feature from amorphous silicates is clearly 
seen, although it often appears to be structured. In IRAS 
17074-1845, 17203-1534, and 20572-^4919 we clearly dis- 
tinguish the presence of single features contributing to 
the overall bump. 

The shift in shape from the smooth ISM-like feature 
(clearly peaked at 9.8 /im) to the structured plateau- 
like feature of the material in the Solar System (peaked 
at 11.3 /im) has been explained as due to grain pro- 
cessing, i.e. coagulation of small (~0.1 /im) grains 
into large (1-2 /im) grains and annealing of amorphous 
silicates into crystallin e magnesium silicates and silica 
(|Bouwman et al.ll2001l ). 

A correlation of the shape and strength of the 10- 
/im silicate emission feature has been observed in Her- 
big AeBe and T Tauri stars and it has been interpreted 
as evidence of grain processing iii the circumstellar disks 
of th ose stars (jvan Boeckel et all 120031 : iPrzvgodda et al.l 
[200l . 

The shape of the feature can be estimated by com- 



paring the continuum-subtracted flux ratio at 11.3 and 
9.8 /im to the peak-to-continuum ratio. We have per- 
formed this comparison for our 0-rich targets. Mixed- 
chemistry sources were not included because of the PAH 
11.2 /im feature, which would affect the flux estimates. 
IRAS 18435-0052 was not included either, because it 
shows the 9.8 /im feature in absorption, which would 
clearly make it an outlier. The plot obtained is shown 
in Figure [HI The data points are sparse in the diagram: 
though a similar trend as in young stars might be present, 
we cannot conclude that this is the case. The absence of 
linear correlation may be explained as due to non-linear 
effects in highly p rocessed dust (11.3 /9.8 flux ratio > 1), 
as pointed out bv lApai et al.l (120051) . O ur results match 
with those found bv lGielen et all (|2008f ). who inspected 
disks around post-AGB stars. The similarity of the re- 
sults leads us to the conclusion that in spite of the ab- 
sence of correlation, the dust in our targets shows a high 
degree of processing, which points to the presence of a 
stable structure where the grains can have time to grow 
and anneal. This is a hint for the presence of circumstel- 
lar/circumbinary disks in our sources. 

As noticed above, IRAS 18435-0052 is the only target 
in which the silicate bump shows an absorption feature 



8 



Cerrigone et al. 




6 8 10 12 



15 20 25 30 
A [/iJ-rn] 



Fig. 6. — As in Figure[3l but with sources classified as mixed for the presence of both PAH and silicate features. As an example, error 
bars are overplotted for IRAS 19336-0400. 



on top of the emission. Since silicates are seen in absorp- 
tion early after the AGB, the presence of both absorption 
and emission features may be interpreted as an earlier 
evolutionary stage for this source. Althoug h this tar- 
get h as been classified as a BII star (Parthsa rathv et al.l 
120001 ) ■ its IRS spectrum indicates major differences to 
the rest of the sample. Almost all of these targets also 
show an emission feature around 18 /im that can be at- 
tributed to amorphous silicates. In Figure [10] we show 



two examples of this feature. 

In such targets as 17381-1616, 17460-3114, 19157-0247, 
18371-3159, 19336-0400, 19590-1249, 20462-1-3416, and 
22495-1-5134 crystalline silicate features around 23.7, 27.6 
and 33.6 /xm are visible. The presence of crystalline sil- 
icates in post-AGB stars has been previously reported 
although different hypotheses exist about the pro duction 
mechanism of these crystals ()Molster et al.|[l999() . 

In Figures [6] and [7] the spectra of those sources with 
both silicate and PAH features are shown. The 10 fj,m 
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silicate bump is typically much stronger than the 6.2 and 
7.7 /im PAH features. The 11.2 fim feature is blended 
with the silicate emission and detected on top of it. Six 
of these stars also show the [Ne II] 12.8 ^m line, which 
proves their ionized status. Four of them also have the 
[Ar II] 7 /im emission line. 

3.3. PAH features 

The presence of PAH features in post-AGB envelopes 
is well-known, although it is not yet clear how these 



molecules are produced ('Tielens"2008D. 

As a first step in the analysis of these features, we 
have classi fied our target s according to the scheme de- 
veloped in iPeeters et al.l ()2002f ). with the exception of 
IRAS 06556-^1623 and 17542-0603 because of their S/N 
ratio in the 5-11 /im range. Almost all of our sources 
have their peak of the 6.2 /im emission feature beyond 
6.235 /im, which puts them in class B. If we also consider 
the 7.7 /im and 8.5 /im features, the stars in our sam- 
ple show a more complicated picture. While the three 



10 



Cerrigone et al. 



lEAS 06556+1623 



0.6 - 




12 14 16 18 20 22 24 
X [/Lira] 



Fig. 8. — Feature around 17—18 /im detected after subtraction of a 
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Fig. 9. — Grain processing (flux at 11.3 over flux at 9.8 fim) as a 
function of the 10-fim silicate emission feature strength (peak flux 
over continuum flux). 

A class sources are also classified as A' and A", class B 
sources spread into a group of intermediate subclasses. 
Tables [H and [3] report the classification for every tar- 
get. The wavelengths calculated do not correspond to 
the peaks of the features, but to their central (median- 
flux) wavelengths, which better accounts for the different 
components within each complex. The 7.7 fim features 
have been classified as A' or B' if their central wavelength 
was respectively smaller or larger than 7.7 /.tm, since the 
complex appears to be mainly due to two features peak- 
i nK around 7.6 a nd 7.8 /xm. 

iPe eters et all (pOOl ) propose that objects in class B 
contain "fresh" PAH molecules, since the post-AGB 
stars and PNe they observed fall in this class. Class 
A objects, which include PDRs and reflection nebulae 
(ISM- like PAHs), are sources with more processed PAHs. 

In the attempt to characterize the PAH populat ion 
in our targets, we have followed iKeller et all (120081 ) in 
their inspection of PAH features in Herbig AeBe stars. 
We fit a spline to the continuum in the 5-14 fim region 
of the spectra, including in the spline the 10 /im bump 
from amorphous silicates, when present. The continuum- 
subtracted spectra are shown in Figure [TlJ We notice 
that despite the spline subtraction, which should remove 
the 10 fim silicate feature, residual emission in the 9- 
11 /xm range remains. As noted above, these features 
may be due to hot crystalline silicates. 

After calculating the central wavelengths as in their 
study, we compare our results to IKeller et all (|2008D . who 
found a correlation between the effective temperature of 
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locally-deflned continuum in two C-rich envelopes. 

the central star and the central wavelength of both the 
11.2 and the 7.7 /im features. Although the astrophysical 
targets are quite different (young and post-AGB stars), 
the properties of the PAHs can be compared to deter- 
mine how these molecules behave in different environ- 
ments. In Figure [121 we show the correlation between 
the central wavelengths of the features at 6.2, 7.7, and 
11.2 fim. A trend with increasing temperature from red 
to blue wavelengths is clear. With the exception of two 
outliers (22495-^5134 in the top and 18379-1707 in the 
bottom of the 7.7 vs 11.2 / zm plot), our data follow the 
trend in I Keller et all (|2008f ) , being located at the high- 
temperature limit, as expected for stars with an ionized 
envelope (Te// > 20000 K). Although the plots show a 
clear correlation, the position occupied by the hot PN 
NGC 7027 (Te// ~ 10^ K) - shown as a dark green dot 
just as an example of hot PN - points out that other 
factors play an important role. 

This complicated dependence is evident when inspect- 
ing the PAH ionization fraction. Since the intensity of 
the features in the 6-9 fj,m range relative to that of the 
11.2 /im one is an order of magnitude higher for ion- 
ized than for neutral PAHs (jGalliano et al.l 120081 ) , in 
Figure [T3l w e plo t I7.7/I11.2 vs I6.2/I11.2 ratios. While 
Keller et ai (120081 ) could not conclude that a correlation 
was present in this diagram, the inclusion of our sam- 
ple allows us to find a remarkable correlation, but no 
clear link between the ionization fraction and Te//. Our 
sample of post- AGB stars contains hotter stars than the 
HAeBe stars in IKeller et all (|2008l ) (Figure [IS and the 
presence of radio continuum) but it has a somewhat lower 
fraction of ionized PAHs. This trend for a decrease of 
io nization fraction in t he post-AGB has been reported 
in iMolster et al.l (|199(]f) for 3 objects. We find this as a 
general trend within our sample. 

Since the ionization fraction does not seem to corre- 
late with the temperature of the central star, we tried 
to link it to the local physical conditions experienced 
by the PAHs, by plotting the I7.7/I11.2 ratio versus the 
ionization factor GoUeT'^'^, where Gq is the radiation 
field in Habing units, Ug the electron density in cm~'^, 
and T the gas temperature in 10^ K. We have calcu- 
lated this parameter by applying the empirical formula 
in Galliano et al. (2008) that links it to I6.2/I11.2 in its 
Lorentzian formulation, since it extends to lower val- 
ues of I6.2/Iii.2- In Figure [131 we show the result for 
400 <GoneTO-5 < 4000 (limits of rehability of the for- 
mula). Although the data are spread in the plot, a weak 
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Fig. 10. — Feature around 17-18 fim detected after subtraction of a locally-defined continuum in two O-rich envelopes. 



TABLE 2 

Central wavelengths of the P AH features and rela tive classification 

ACCORDING TO IPeETERS ET AlI (j2002l '). 



Target 
IRAS ID 


-^6.2 
fim 








At. 7 






As. 5 






All. 2 
fj.m 






Class 


01005+7910 


6.231 


± 





005 


7 


58 







08 


8.59 







03 


11.21 


± 





07 


AA 


A" 


09470-4617 


6.227 


± 





005 


7 


66 


± 





03 


8.60 







06 


11.29 


± 





01 


AA 


A" 


18442-1144 


6.238 


± 





001 


7 


64 


± 





01 


8.60 


± 





07 


11.205 


± 





004 


BA 


A" 


19200+3457 


6.24 


± 





04 


7 


61 


± 





07 


8.58 







02 


11.250 


± 





006 


BA 


A" 


18371-3159 


6.28 


± 





01 


7 


68 


± 





09 


8.61 


± 





09 


11.2 


± 





1 


BA 


A" 


18379-1707 


6.266 


± 





004 


7 


62 


± 





04 


8.57 


± 





06 


11.03 


± 





06 


BA 


A" 


19306+1407 


6.240 


± 





003 


7 


74 


± 





04 


8.6 


± 





1 


11.30 


± 





07 


BB 


B" 


19336-0400 


6.25 


± 





01 


7 


61 


± 





02 


8.70 


± 





02 


11.30 


± 





03 


BA 


B" 


19590-1249 


6.22 


± 





06 


7 


58 


± 





02 


8.62 


± 





03 


11.21 


± 





03 


AA 


A" 


20462+3416 


6.25 


± 





01 


7 


66 


± 





06 


8.58 


± 





04 


11.26 


± 





03 


BA 


A" 


21289+5815 


6.251 


± 





001 


7 


74 


± 





03 


8.64 


± 





08 


11.25 


± 





07 


BB 


B" 


22023+5249 


6.23 


± 





01 


7 


69 


± 





03 


8.55 


± 





09 


11.23 


± 





07 


BA 


A" 


22495+5134 


6.25 


± 





01 


7 


58 


± 





02 


8.5 


± 





1 


11.45 


± 





05 


BA 


A" 



TABLE 3 
Intensity of the PAH features. 



Target 


l6.2 


I7.7 


Is 


5 






111. 2 


IRAS ID 


10^15 


10-15 W m-2 


10-15 W 


m 


-2 


10-15 W m-2 


01005+7910 


24 ± 1 


38 + 2 


9 


± 


1 




29 ± 3 


09470-4617 


11.5 ± 0.6 


19 + 1 


2.6 







5 


14 ± 1 


18442-1144 


4.34 ± 0.05 


8.2 ± 0.3 


2.2 


± 





2 


2.41 ± 0.04 


19200+3457 


1.77 ± 0.05 


3.3 ± 0.1 


0.64 


± 





06 


1.10 ± 0.03 


18371-3159 


0.12 ± 0.02 


0.25 ± 0.05 


0.04 







01 


0.13 ± 0.05 


18379-1707 


4.0 ± 0.2 


5.6 ± 0.4 


0.5 


± 





1 


1.1 ± 0.3 


19306+1407 


3.98 ± 0.03 


8.6 ± 0.4 


0.2 


± 





1 


4.9 ± 0.8 


19336-0400 


0.17 ± 0.03 


0.74 ± 0.05 


0.06 


± 





01 


0.33 ± 0.04 


19590-1249 


0.17 ± 0.05 


0.37 ± 0.05 


0.03 


± 





01 


0.18 ± 0.06 


20462+3416 


0.49 ± 0.02 


1.5 ± 0.1 


0.14 


± 





03 


0.28 ± 0.07 


21289+5815 


1.20 ± 0.02 


2.2 ± 0.2 


0.18 


± 





05 


0.35 ± 0.07 


22023+5249 


2.4 ± 0.1 


4.7 ± 0.2 


0.6 


± 





1 


1.0 ± 0.2 


22495+5134 


0.31 ± 0.04 


0.52 ± 0.04 


0.13 


± 





03 


0.6 ± 0.06 



correlation can be identified, indicating that the ioniza- 
tion fraction increases with the ionization factor. Once 
again, no clear correlation with Tg// is observed, indi- 
cating that the local conditions in the envelope play a 
major role. 

The results of our analysis of the PAH features are 
somewhat conflicting. On one side we see that the central 
wavelengths of the features are shifted toward the blue, 
as expected in photo-processed molecules; on the other 
side, we find a low ionization fraction, which implies low 



processing. 

We conclude that a higher effective temperature im- 
plies a harder field and stronger processing in terms of 
shift of the central wavelength toward the blue, but not 
necessarily a higher ionization fraction, this being depen- 
dent on the available amount of ionizing photons, thefore 
on the distance to the central star. Our plots indicate 
that the PAH molecules in our stars are located in the 
outflows, where they are subject to a more diluted field. 
Although the probability of interaction between photons 
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Fig. 11. — Features in the 5—13 fim range after subtraction of a spline continuum that includes the silicate bump. The first five digits 
of the IRAS ID are listed next to each plot; in the bottom right an example of a spline fit is also shown, with anchor points as diamonds. 
The vertical dotted lines indicate the typical peak wavelengths for PAH features. 



and molecules far from the star is smaller, the hard pho- 
tons from the central star - when interacting - determine 
a major damage to the molecules, resulting in the shift in 
central wavelength. In this scenario, the closer molecules 
must have been previously destroyed by the radiation 
field. Since PAH features with red central wavelengths 
are found in post-AGB stars (RAFGL 2688 i s one of the 
two class C objects in IPeeters et al.l (j2002D ). this dis- 
ruption should occur during the few 10'^ yr between the 
early post-AGB and the development of an ionized shell. 
High-angular resolution observations in the PAH emis- 
sion ranges are necessary to support such an interpreta- 
tion, which confirms the idea of C-bearing molecules be- 
ing located in the outflows of mixed-chemistry nebulae. 
Finally, it must be noticed that the wavelength shifts ob- 
served can be determined not only by photo-processing: 



nitrogenation can in fact be another explanation (jTielensI 
I20n8[) . 

3.3.1. Abundance estimate 

As shown in lAllamandola et aTl (|1989D , the mass that 
goes into PAHs can be linked to the ratio of the flux 
emitted in the PAH features to the continuum flux in 
the far infrared range (FIR), nominally between 40 and 
500 /im. We have therefore taken advantage of our SED 
models (see § 2]) to calculate the FIR flux by integrating 
our curves in the mentioned interval. 

The fra ction of carbon trapped in PAHs can be esti 
mated as (|Allamandola et al 



pped n 

MM) 



fc 



1.6 X IQ- 



-18 



(1) 
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Fig. 12. — Corr elation between the central wavelengths of the 6.2 and 11.2 features with that of the 7.7 complex. The data from 
IKeller et all 1)20081 ) are shown with a color proportional to the stellar temperature ranging from red (cold) to blue (hot). Our data are 
shown in yellow (mixed chemistry) and blue marine (C-rich). NGC 7027 is shown in green. 
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Fig. 13. — (a): The intensity ratio of the 7.7 to 11.2 /im feature versus the 6.2 to 11.2 ratio. A clear correlation is observed, which 
points to the dependence of these ratios on the ionization fraction. No clear correlation with the temperature of the central star is found, 
(b): The I7.7/I11.2 ratio, pr opo rtional to the ionization fraction, plotted versus the ionization factor, which accounts for the local physical 
conditions. Colors as in Fig [12] 



where auv is the average UV absorption cross section of 
PAHs per C atom and fir is the PAH to continuum flux 
ratio. We assume 7 x 10~^^ cm~^ as a typical value for 
([Tielens 2005). 
The fraction of C in PAHs can then be estimated as: 



fc = 0.23 



fit 



1 - /« 



(2) 



Table [3] lists the results of the calculation, along with 
the estimates of the number abundance of PAH molecules 
to H nuclei, assuming that the average PAH molecule 
has 50 C atoms. The errors in the table have been cal- 
culated by standard error propagation, including a 15% 
uncertainty in the model estimation of the underlying 
continuum. 

The results show a trend for larger values of fc and 
PAH abundance in C-rich stars. In these sources, the 
values of C fraction range from 1% to 43% (3% is a typical 
value in the ISM), while those of the PAH abundance are 



between 7x10"^ and 336x10"^. With the remarkable 
exception of IRAS 21289-^5815, stars with both PAH and 
silicate features show lower values: 0.1-1% for fc and 
0.5-8 X 10~* for the PAH abundance. We notice here that 
since in our modeling we lack observations in the sub- mm 
and mm ranges, our estimation of the FIR flux and the 
values derived from it should be taken with caution. 

4. MODELING THE SPECTRAL ENERGY DISTRIBUTION 

To model our targets, we have collected data from 
the 2MASS (jSkrutskie et all I20060 and IRAS archives. 
We have dere d dened the 2MASS magnitudes following 
iSchlegel et al.l (jl998). The intrin sic B -V color indexes 
have been estimated from Lang (|1999( ) on the basis of 
the spectral classifications of the targets, then the mag- 
nitude correction was calculated with the derived color 



We have used the DUSTY code (jlvezic et al.lll999D to 
solve the problem of radiation transfer in the dust en- 
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TABLE 4 

Fraction of C trapped in PAHs fc (column 2) and 
PAH abundances (column 4). 



Target 










IRAS ID 


10^2 


10-2 


io-« 


10-* 


C-rich 










01005+7910 


6 


2 


51 


19 


09470-4617 


43 


37 


336 


293 


18442-1144 


1.0 


0.3 


8 


2 


19200-f3457 


2.7 


0.9 


21 


7 


Mixed 










18371-3159 


0.07 


0.02 


0.55 


0.18 


18379-1707 


1.1 


0.4 


8.7 


2.8 


19306-fl407 


0.48 


0.15 


3.7 


1.1 


19336-0400 


0.12 


0.04 


0.9 


0.3 


19590-1249 


0.10 


0.03 


0.7 


0.2 


20462+3416 


0.12 


0.03 


0.9 


0.3 


21289+5815 


5.2 


1.9 


41 


15 


22023+5249 


0.5 


0.1 


3.8 


1.2 


22495+5134 


0.12 


0.04 


1.0 


0.3 


Note. — Columns 


3 and 5 


list the errors 


associated to 



fc and ^pAH /^H calculated including a 15% uncertainty 
in the estimate of the underlying continuum. 

velopes of our targets. 

Our typical first cycle of DUSTY had the following 
input: 

- central source: blackbody curve with temperature 
matching the spectral type of the target; 

- chemistry of the envelope: either 100% silicates or 
100% amorphous carbon, depending on IRS fea- 
tures; 

- grain size distribution: power law oc a~^'^, a-min = 
0.005, araax = 0.26, a being the grain radius in /im; 

- density distribution: power law oc i?^^, R being 
the shell radius; 



shell relative thickness: Ro 



1000 Rr. 



- dust temperature at Rm'- adapted to each SED; 

- optical depth at 60 /im: calculated for each star 
from IRAS data; 

We calculated the opacity at 60 /xm for each object as 



Tu 2.03 X 10 



10 



(3) 



where is the dust temperature, the angular radius 
of the nebula in arcsec, Fi, the flux density at the fre- 
quency v in erg cm"^ s"^ Hz"^, Bi,{Td) the Planck func- 
tion at the tempera ture in erg cm"^ s"^ Hz"^ sr"^ 
()Gathier et al.lll98l . This gives us a first estimate of the 
opacity, based on the 5" upper limit to the envelope sizes 
derived from our IRAC images. This value of opacity was 
used in the first cycle of DUSTY and then changed in the 
following cycles to better match the data. 

In several targets we have found that it is not possi- 
ble to model the SED with only one dust component. 
DUSTY does not allow the user to include in the model 
several dust shells at the same time, but it is possible to 
use the output of a run as an input to a second run, thus 



approximately reproducing a multi-shell structure. We 
have successfully applied this technique to the sources 
that clearly showed a near-IR excess of radiation, point- 
ing out the presence of hot dust. A summary of the in- 
put parameters used in our modeling is given in Table El 

One of the issues encountered in fitting these data in- 
volves the actual shape of the 10 /im bump. Our models, 
though often reproducing the overall shape of the SED 
quite well, fail in reproducing the correct shape of the 10 
/tm bump especially in its longer wavelength part, where 
it overlaps with the underlying rising continuum. The 
shape of the feature depends on several factors, including 
dust temperature, grain size, and composition. Geomet- 
rical effects can also play a role, if different dust species 
are segregated in different regions of the envelopes. 

Another issue concerned the fit of the near-IR emis- 
sion. In several instances DUSTY underestimates the 
continuum in this range. This can point out the pres- 
ence of additional hot dust, but another possibility is 
the presence of a free-free contribution from the ionized 
envelope, which is neglected in DUSTY. Our radio obser- 
vations prove to be very useful then, since by modeling 
the radio emission we can check if the hypothesis of a 
free-free contribution in the near-IR is valid or not. 

We have modeled the radio continuum spectra assum- 
ing the central star is surrounded by a shell of ionized 
gas with constant density. The usual approximation of 
the Gaunt factor at radio wavelengths typically results 
into a quick cutoff of the free-free emission in the far-IR, 
where, compared to the dust emission, such contribu- 
tion is indeed negligible. Nevertheless, in the near-IR, 
where the central star and dust fluxes are at a minimum, 
the free-free contribution to the observed flux may not 
be much smaller than the contributions from the other 
emitting components. For example, in the evolved PN 
IC 4406 a near-IR excess, previously supposed to in- 
dicate a hot-dust com ponent, has been expla ined as a 
free-free contribution (jCerrigone et al.l l2008bD . There- 
fore, we decided t o esti mate the Gaunt factor following 
iKarzas fc Latteij (11961) for a 10** K electron tempera- 
ture, extending the calculation of the emission curve from 
the cm ra nge to 0.8 /im. F o r detai ls about the radio mod- 
eling, see iCerrigone et all (|2008a( ). 

Where radio data were available, we summed up the 
free- free and dust model contributions. This often results 
into a substantial correction to the dust-only model, not 
only in the cm range, as expected, but also in the near-IR. 
In Figures IHl and 1T51 we show examples of the performed 
DUSTY and radio fits for about half of our sample. 

During the various trials performed to match the data, 
we have found that in some sources it is necessary to 
add an extra-component that does not match chemi- 
cally with the classification of the target. Specifically, 
in half of the 0-rich targets (IRAS 17364-1238, 17542- 
0603, 18435-0052, 19157-0247, 20572-^4919) we had to 
include an am-C component to match the observational 
points. Similarly, in two of the C-rich sources (IRAS 
06556+1623 and 19200-K3457) we had to include a sil- 
icate component. It must be considered that this may 
be due to limits in the modeling code, in particular the 
current knowledge of the opacity coefficients in the mid- 
IR, and does not necessarily represent a mixed-chemistry 
nature of the targets. IRAS 17423-1755 has been mod- 
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TABLE 5 

Parameters used in the DUSTY models. 
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" T* is the blackbody temperature used for the central source. 

^ Sil-Ow, Sil-Oc, Sil-DL, and am-C indicate respectively Ossenk opf warm silicates, Ossenkopf 
cold silicates, Draine-Lee silicates, and amorphous carlion. See llvezic et al.l l|1999l ) for refer- 



Minimum and maximum grain radii used in modeling. 



eled with amorphous carbon grains, ahhough it does not 
show any PAH features. This target is reported to show 
C2H2 absorption feature (jGauba fc Parthasarath-vl2004f ). 
which indicates it has a C-rich chemistry. Such a chem- 
istry also matches with its featureless mid-IR spectrum. 
In our attempts to model it with silicates, we have found 
that if silicates were present they would determine a dis- 
tinct 10-/im emission feature, because of their high tem- 
perature (peak in IRAC bands). In two targets showing 
both PAH and silicate features (IRAS 18379-1707 and 



21289-1-5815) a silicate-only chemistry is sufficient to fit 
the continuum spectrum. 

IRAS 17423 -1755 is our only target with mm measure- 
ments (Huggin s et al.l 120041. We have assumed as size 
of the dust region the angular size of its CO shell as re- 
ported in Huggins e t al. (2004). ~2", with a distance of 
5.8 kpc and luminosity 12600 Lq. When using in DUSTY 
the opacity calculated at 60 /im from the IRAS data, we 
find that the model highly overestimates the far-IR emis- 
sion. Changes in the opacity input value do not improve 
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Fig. 14. — Examples of the performed SED modeling. The IRS spectra are overplotted in red. Radio data are shown as solid circles, 
IRAS data as squares, IRAC as triangles, and 2MASS as diamonds. On the right, the dashed lines are the sole DUSTY output, whereas 
the solid line is the sum of the DUSTY and the free-free model. 



the model, nor do different density distributions (we tried 
r~^ and r~'^). A good match can be obtained if a maxi- 
mum grain size of 100 /im is introduced fFigure fTB)) . This 
is remarkable, since hints for such large grains have been 
found so fa r only in a handful of stars, such as the Red 
Rectangle (jJura et al.ll200l]) . It is interesting to notice 
that in another target, IRAS 20462-1-3416, our point at 
1 cm seems not to follow the free-free curve, being well 
above it (see Figure [Ti|) . If we combine this to the shape 
given by the IRAS long-wavelength data points, we can 
speculate that this target might have an excess of radia- 
tion in the sub-mm as well. 

We could not classify IRAS 18367-1233 because the 
IRAS coordinates turned out to be off-set respect to the 
actual position of the target, as we can see in our IRAC 
images. Since the IRAC imaging was not available at the 
time of the IRS observations, the target resulted to be 
out of the IRS slits. 

4.1. Circumstellar mass calculation 

Modeling our sources allows us to estimate the mass 
of the en velopes of our targets. We follow here the pro- 
cedure in lSarkar fc Sahail (|2006f) . The mass is estimated 
as 

M = A7TRt^Yinoo/kioo)S (4) 



In Eq. [l]i?i„ is the inner radius of the envelope, Y is the 
outer-to-inner radius ratio, tioo and kiQo are the optical 
depth and the absorption coefficient at 100 fim, and 6 is 
the dust-to-gas mass ratio. No distance information is 
available for most of our targets, therefore we assume as 
Rin a typical value of 10^^ cm. Y is the value given as in- 
put in DUSTY and tioo is given by the code as an output. 
For S, whose value depends on the evolutionary stage, 
we consider a typical value of 200. fcioo depends on the 
chemistry of t he envelope: we use 92 g cm~^ for amor- 
phous carbon ( Cerrigone et al. 2008b) and 34 g cm for 
silicates (Sa rkar fc Sahai 2006). We have used the values 
for silicates not only in the targets that we have classi- 
fied as 0-rich, but also in mixed-chemistry environments, 
neglecting any contribution fro m C-b e aring molecules. 
As already noticed by iSarkar fc Sahail (j2006l ) , the code 
seems to underestimate the amount of cold dust in the 
envelope or anyhow not to be very sensitive to the contri- 
bution from such grains. If we also consider that we lack 
measurements in the sub-mm/mm range - which would 
help constrain the cold-dust emission - the values listed 
in Tablets] are to be regarded as lower limits to the cir- 
cumstellar mass. 

5. SOURCES WITH UNCLASSIFIED SPECTRA 
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Fig. 16. — SED and DUSTY model o f IRA S 17423-1755. The IRS spectrum is overplotted in red. Radio data are shown as solid circles 
(1 and 3 mm data from lHuggins et al.l II2004I) ). IRAS data as squares, IRAC as triangles, and 2MASS as diamonds. 
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TABLE 6 

Estimates of the 
envelope masses. 
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For IRAS 18040-1457, 18070-2346, and 19399-1-2312 
the IRAS measurements must evidently include emission 
from nearby dust, since in our IRAC images we clearly 
detect a diffuse nebulosity around these targets (see as an 
example Figure I17p . These sources have quite different 
IRS spectra, compared to the other targets (Figure [TH]) . 
Since IRAS data are likely contaminated by ambient 
dust, the location of these targets in the IRAS color- 
color diagram is probably wrong, therefore their clas- 
sification as post-AGB candidates appears dubious. In 
IRAS 18070-2346, the 18 /.tm silicate feature is detected, 
although it might be due to the interstellar medium or 
diffuse nebulosity. IRAS 19399+2312 shows an emission 
feature around 11.2 /im that can be attributed to neutral 
PAH molecules. 

6. CONCLUSIONS 

The number density of post-AGB/proto-PN is esti- 
mated to be around 0.4 kpc^^. As a comparison, the 
density of AGB stars is 15 kpc~^ and of M ain Sequence 
stars is -2x10*^ kpc'^ (| Woods et all 120051 ). Knowledge 
about this phase is mostly derived from a handful of ob- 
jects, such as the C-rich source RAFGL 618 and the 
O-rich object OH231.8-f 4.2. Because of the paucity of 
known TO, the identification and study of new TO is 
very important for testing models of stellar evolution. 

Among TO candidates, we have selected hot targets (B 
spectral type), therefore searching for recently- ionized or 
close-to-ionize envelopes. 2MASS and IRAS data were 
combined with our IRAC and radio observations to con- 




FlG. 17.— IRAC image of IRAS 18040-1457 at 8 ^tm in logarith- 
mic scale. The field is about 3'. 6 X 3'. 6 and the units are MJy/sr. 
North is up and East is left. 



struct the SEDs of the sample stars. DUSTY modeling 
shows the presence of more than one emitting component 
in the dust envelopes and indicates that dust tempera- 
tures range from ^1000 K to ^100 K. The modeling 
also gives hints about the presence of large (> 1 /xm) 
grains, although sub- mm and mm observations are nec- 
essary to constrain the SEDs in these ranges, where such 
grains would give a major contribution. As an example, 
in IRAS 17423-1755 it is necessary, to account for its mm 
emission, to have three dust components (at 1150, 132, 
and 80 K) with grain size up to 100 /im (see Fig [TO]) . 

The IRS spectra we obtained allow us to classify the 
TO candidates as C-rich, O-rich, or mixed-chemistry. We 
have classified as mixed-chemistry objects those stars 
with both a 10 /im silicate bump and PAH features 
in the 6-11 /im range. A possible explanation for the 
mixed chemistry involves the presence of a circumbinary 
disk/torus, where 0-bearing molecules would be pre- 
served from the third dredge up, while elsewhere in the 
outflow features from C-bearing molecules would arise. 
This picture appears to be confirmed by our analysis of 
the PAH features, which indicates that PAH molecules 
in our targets are not located close to the central star. 
Since these objects are typically compact (^--^2"), sub- 
arcsec imaging may evidence if different dust species have 
different spatial distributions. 

In most of our targets with O-rich dust, the 10 fxin 
silicate emission feature often appears with a structured 
shape and/or a peak shift to larger wavelengths. Our 
analysis of the shape and strength of this feature indi- 
cates a high degree of dust processing. Such process- 
ing, which typically consists in dust growth and crys- 
tallization, indicates again the presence of large dust 
grains and/or of crystalline silicates and more interest- 
ingly is a further hint for the presence of a circum- 
stellar/circumbinary disk: a stable structure where the 
grains can have the time to grow and crystallize. 



PAH and silicates in post-AGB stars 



19 



0.6 r 



0.5 



IRAS 18040-1457 




0.0 , 



IRAS 18070-2346 



2.0 



1.5 



1.0 - 



0.5 



0.0 ; 



0.3 



0.2 - 



0.1 - 



0.0 [ 




10 



15 



20 25 
A [/d.in] 



30 



35 



Fig. 18. — Spectra of sources with dubious classification. Error bars are overplotted. 



Whereas the expected fraction of mixed-chemistry en- 
velopes is less than 10%, we find that about 40% of the 
stars in our sample show both PAH and silicate features. 
Our sample has been selected on the basis of B optical 
spectral type and far-IR excess. These two conditions are 
indeed typical of hot post-AGB stars. Nevertheless, se- 
lections based on infrared thermal emission (IRAS fluxes, 
for instance) introduce a bias towards far-IR bright ob- 
jects. Among post-AGB stars, objects hosting a disk 
or torus are stronger IR emitters. These circumstellar 
(or circumbinary) structures are highly stable and store 
the dust around the central star/s on a much longer 
time-scale than in stellar envelopes without a disk/torus 
(jvan Winck el 2003) . As these envelopes are brighter and 
live longer, selections based on far-IR emission are bi- 
ased towards stars surrounded by a disk or torus. The 
higher percentage of mixed-chemistry envelopes in our 
sample can therefore be connected to a higher percentage 
of disks (or tori) . This confirms the general picture of a 



disk/torus as a reservoir preserving 0-bearing molecules 
from the third dredge up. 

If we refer to the Red Rectangle as an example, mixed 
chemistry and large grains are linked together by the 
presence of a disk, where the particles can survive enough 
orbits to grow by coagulation to sizes as large as 1 mm 
(jJura et al.lll997l ) . With only a few examples of this phe- 
nomenon, several hints for it found in this project provide 
us with candidate targets to test the current understand- 
ing of dust processing in PN envelopes. 
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